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Chapter 1 

Gaussian Behavior in Generalized Zeckendorf 

Decompositions 

Steven J. Miller and Yinghui Wang 



Abstract A beautiful theorem of Zeckendorf states that every integer can be writ- 
ten uniquely as a sum of non-consecutive Fibonacci numbers {Fn}'^^^; Lekkerk- 
erker proved that the average number of summands for integers in [F„,F„+i) is 
n/{(p^ + 1), with (p the golden mean. Interestingly, the higher moments seem to 
have been ignored. We discuss the proof that the distribution of the number of 



(~| , summands converges to a Gaussian as « -^ oo, and comment on generalizations 



to related decompositions. For example, every integer can be written uniquely as 
a sum of the ±F„'s, such that every two terms of the same (opposite) sign differ 
in index by at least 4 (3). The distribution of the numbers of positive and negative 
summands converges to a bivariate normal with computable, negative correlation, 

namely -(21 -2^)/(29 + 2^)w -0.551058. 



Keywords: Fibonacci numbers, Zeckendorf's Theorem, Lekkerkerker's theorem, 
r — ] generating functions, partial fraction expansion, central limit type theorems, far- 

^SJ ' difference representations. 

(3 ■ MSC 2010: 1 1B39 (primary) 65Q30, 60B10 (secondary). 



1.1 Introduction 
c^ ' 1.1.1 History 



The Fibonacci numbers have fascinated professional mathematicians and ama- 
teurs for centuries. They have a wealth of properties and interesting relationships; 
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see for example IKosl . This article is concerned with how Fibonacci numbers arise 
in decompositions. A beautiful theorem of Zeckendorf QZel states that every posi- 
tive integer may be written uniquely as a sum of non-adjacent Fibonacci numbers; of 
course, to ensure that the decomposition is unique we need to use the normalization 
Fi = 1, F2 = 2, Fi, ~3, F4~5 and in general /Vi+i = F„ +F„_i. 

The standard proof is by induction, and is also constructive. Simply take the 
largest Fibonacci number that is at most our number x; say this is F„. If F„ +F„_i < x 
then we could take F„^ 1 in our decomposition, and then replace Fn and F„^ 1 with 
Fn+i by the recurrence relation, which contradicts the maximality of F„. Thus 
F„ + f„_ 1 > X and F„_ 1 cannot be used. By induction, we can write x — F„ as a 
sum of non-adjacent Fibonacci numbers; further, since x — F„ < F„^ \ clearly F„^i is 
not in our decomposition. Adding F„ to the decomposition of x — F„ yields our de- 
sired expansion. To prove uniqueness, we again proceed by induction. Arranging the 

Fibonacci summands in decreasing order, if x = F„, H \-F„i, = F„^^ h/^,np we 

must have F„, = F,„, ; if not, the two sums cannot be equally The claim now follows 
by induction as x — F,^ < x. 

There are many questions one may ask about the Zeckendorf decomposition. 
The first is to understand the average number of summands needed. Lekkerkerker 
BLekl proved that for x € [Fn,Fn+i), as « -> 00 the average number of summands 

needed is «/((p^ + l), with (p = ' ^ the golden mean. The proof involves continued 
fractions, and is unfortunately limited to the mean (though it can be generalized to 
other related decompositions; see IBCCSW Day , Ha. HojIKellLenl for some of the 
history and results along these hnes). 

In this chapter we discuss the combinatorial vantage of Kologlu, Kopp, Miller 
and Wang that allows us to determine, not just the mean, but all the moments and 
thus the limiting distribution. They prove the fluctuations of the number of sum- 
mands needed forx G [F„,F„+i) about the mean converges to a Gaussian as n -^ °o. 
The proof follows from writing down an explicit formula for the number of x where 
there are exactly k summands in the Zeckendorf decomposition, and then using 
Stirling's formula to analyze the resulting binomial coefficients. These results were 
generalized in Miller- Wang [|MW| to other decompositions in terms of elements 
satisfying special recurrence properties. While the combinatorial approach is still 
the foundation of the analysis, the resulting expressions are too involved and can- 
not be attacked as directly as the Stirling approach in the initial case. Instead, the 
proofs are completed by an analysis of the resulting generating functions, and then 
differentiating identities to determine the moments. 

We describe the main results in i)1.1.2| and then sketch the proofs in the Fibonacci 
case in m.2\ The key idea is recasting the problem from number theory to combina- 
torics, specifically to the number of integer partitions satisfying certain constraints. 
We describe this idea in detail, and outline the resulting algebra (complete details 
are available in QKKMWI ). Unfortunately, while the combinatorial framework ap- 
pUes to generalizations of Zeckendorf decompositions, the proof in general is sig- 



' If F„,| > F„| , then the largest the n-sum can be is F,^ + F„j -2 + ^ni -4 H 1-5, where 5 is either 

1 or 2. Adding 1 or 2 to this and using the recurrence relation gives F„^^i, and thus the expressions 
are unequal. 
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nificantly harder. This is due to the fact that for the Fibonacci numbers, we have 
an exphcit formulas for the probability a number in [F„,F„^i) has exactly ^+1 
summands in its Zeckendorf decomposition; it is just (" ^, ) /fn-i- All our results 
follow from a careful analysis of the behavior as n — > oo, which can be accom- 
phshed with Stirling's formula. In the general case, the resulting expressions are not 
as tractable, and the combinatorial approach must be supplemented. We sketch the 
main ingredients in the analysis in m.3\ and refer the reader to liMW] for full proofs. 
We end in ij 1.41 with some open questions. 



1.1.2 Main Results 



Before stating the main results, we first set some notation. The sequences defined 
below are the generalizations of the Fibonacci numbers that can be handled using 
the techniques of BMWI . 

Definition 1.1.1 A sequence {^n}~=i of positive integers is a Positive Linear Re- 
currence Sequence (PLRS) if the following properties hold: 

1. Recurrence relation: There are non-negative integers L^c\,. . . ,cl such that 

H„+i = ciH„-\ ^clH„+\^l, (1-1) 

with L, c\ and ci positive. 

2. Initial conditions: Hi = I, and for 1 <n<Lwe have 

H„+x = ci//„+C2//„_i + ---+c„//i + l. (1.2) 

A decomposition YULi <^iHm+\-i of a positive integerN (and the sequence {fl(}^jj 
is legal if ai > 0, the other ai > 0, and one of the following two conditions holds: 

• Condition 1: We have m < L and ai = Cifor 1 < / < tn. 

• Condition 2: There exists s € {0, ■ ■ ■ ,L} such that 

a\ = ci, (32 = C2, ••• , fl.s-i = c.s-i anda.v <c,, (1.3) 

a.5+1 , . . . , a^+c = 0/or some i>Q, and {^(}"li* (with bi = ag+t+i) is legal. 

IfYJi'Li'^iHm+i-i is a legal decomposition ofN, we define the number of sum- 
mands (of this decomposition ofN) to be ai -\ \- a„,. 

Informally, a legal decomposition is one where we cannot use the recurrence 
relation to replace a linear combination of summands with another summand, and 
the coefficient of each summand is appropriately bounded. For example, if //„+i = 
2H„ + 3i/„_ 1 + //„_2, then H5 + 2H4 + 3//3 + Hi is legal, while H^ + 2H4 + 3Hj, + H2 
is not (we can replace 2//4 + 37/3+7/2 with 7/5), nor is 77/5+27/2 (as the coefficient 
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of 7/5 is too large). Note the Fibonacci numbers are just the special case of L = 2 
and ci = C2 = 1. 

The following probabilistic language is convenient for stating some of our main 
results. 

Definition 1.1.2 (Associated Probability Space to a Positive Linear Recurrence Sequence) 

Let {H„} be a Positive Linear Recurrence Sequence. For each n, consider the dis- 
crete outcome space 

H„+i-l} (1.4) 



fl C i2„; (1.5) 



in other words, each of the //„+i — //„ numbers is weighted equally. We define the 
random variable K„ by setting K„{(o) equal to the number of summands of (O € i2„ 
in its legal decomposition. Implicit in this definition is that each integer has a unique 
legal decomposition; we prove this in Theorem\l\ and thus K„ is well-defined. 
We denote the cardinality ofQ,, by 

An = H„+i-H„, (1.6) 

and we set /?„<■ equal to the number of elements in [//„,//„_|_i) whose generalized 
Zeckendorf decomposition has exactly k summands; thus 

p„^k = A„-Pmh{K„=k). (1.7) 

Before stating the main results, we first examine a special case which sug- 
gests why the limiting behavior is Gaussian. Consider the PLRS given by L = 1, 
H„+i = IQH„. Thus our PLRS is just the geometric series 1,10, 100, . . . , and a legal 
decomposition of A^ is just its decimal expansion. Clearly every positive integer has 
a unique legal decomposition. Further, the distribution of the number of summands 
converges to a Gaussian by the Central Limit Theorem, as we essentially have the 
sum of « — 1 independent, identically distributed discrete uniform random variables. 

To see this, write N ~ ailQ" -\ h «,,+ 1 1 . We are interested in the large « behavior 

offliH hfl„+i as we vary over x in [ 1 0" , 1 0"^ ' ) . For large n the contribution of a 1 

is immaterial, and the remaining fl,'s can be understood by considering the sum of n 
independent, identically distributed discrete uniform random variables on {0, . . . , 9} 
(which have mean 4.5 and standard deviation approximately 2.87). Denoting these 
random variables by A,-, by the Central Limit Theorem A2 H h A„+i converges to 



being normally distributed with mean 4.5/1 and standard deviation y/33n/4. 

The first result is that a PLRS leads to a unique, generalized Zeckendorf decom- 
position. Part (a) has also been recently studied by Hamlin IHal . 
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Theorem 1 (Generalized Zeckendorf 's Theorem for PLRS IMWl ). Let {//„}~=i 
be a Positive Linear Recurrence Sequence. Then 

(a) There is a unique legal decomposition for each integer N >0. 

(b) There is a bijection between the set S^„ of integers in [i/„,//„_|_i) and the set 
^„ of legal decompositions Y!i=i fl/^n+i-j- 

In addition to being of interest in its own right, the Generahzed Lekkerkerker 
Theorem below is needed as input to prove the Gausian behavior of the number of 
summands. 

Theorem 2 (Generalized Lekkerkerker 's Theorem for PLRS lIMWl ). Let {//„}"= 
be a Positive Linear Recurrence Sequence, let K„ be the random variable of Def- 
inition \LL2\ and denote its mean by jj.,,. Then there exist constants C > 0, d and 
7i € (0, 1) depending only on L and the c,'i in the recurrence relation of the H„ 's 
such that 

lln = Cn+d + o{fi). (L8) 

The main resuh in the subject is 

Theorem 3 (Gaussian Behavior for PLRS OMWI ). Let {//„}"= i be a Positive Lin- 
ear Recurrence Sequence and let K^ be the random variable of Deftnition \LL2] As 
n ^ oa, the distribution ofK„ converges to a Gaussian. 

Their method generaHzes to a multitude of other problems, and allows us to prove 
Gaussian behavior in many other situations. We state one particularly interesting 
situation. 

Definition 1. We call a sum of the ±f„'s a. far-difference representation if every 
two terms of the same sign differ in index by at least 4, and every two terms of 
opposite sign differ in index by at least 3. 



Albert [Al | proved the analogue of Zeckendorf s Theorem for the far-difference 
representation. It is convenient to set 

^ ^ f Eo<n-4K« ^n-4/ = Fn + Fn-4 + Fn-% + ' ' ' if H > 

1 otherwise. 

Theorem 4 (Generalized Zeckendorf 's Theorem for Far-Difference Represen- 
tations UAH ). Every integer has a unique far-difference representation. For each 
N € {S„-i = F„ — 5„_3 — l,S„], the first term in its far-difference representation is 
F„, and the unique far-difference representation ofO is the empty representation. 

The far-difference representations have both positive and negative summands, 
which opens up the fascinating question of how the number of each are related. It 
turns out they are correlated, though the total number of summands and the differ- 
ence in the number of positive and negative summands are not. Specifically, 
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Theorem 5 (Generalized Lekkerkerker's Theorem and Gaussian Behavior for 
Far- Difference Representations HMWI ). Let J(fn and ^„ be the corresponding ran- 
dom variables denoting the number of positive summands and the number of neg- 
ative summands in the far-difference representation for integers in (5„_i,5„]. As n 
goes to infinity, the expected value of ,J^n, denoted by ¥\Jif„], is j^n + Iq 

and ^^ = |- greater than E[^]; the variance of both is of size iqoo »■" 
the joint density of J€„ and ^, is a bivariate Gaussian with negative correlation 
— 175 — ~ ~ 29+2a> ~ —0.55 1; and J^„ + ^„ and J(f„ — ^„ are independent. 
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1.2 The Fibonacci Case 

As the framework of the Fibonacci case is the basis for the general case, and as the 
Fibonacci case can be handled elementarily, we describe it in some detail. The key 
idea is to change our perspective and apply the solution to a standard problem in 
combinatorics, the stars and bars problem. It is also known as the cookie problem in 
some circles (see for instance Chapter 1 of IMT-BI '). or Waring's problem with first 
powers (see IINal ). 

Lemma 1.2.1 The number of ways to partition C identical objects into P distinct 
sets is (*^p^7'). 

Here it is very important that the objects being partitioned are indistinguishable; 
all that matters is how many objects are given to a specified set, not which objects 
are given. 

Proof. Imagine instead C + P— 1 objects. There are { p_^ ) ways to choose P— 1 
of the C. Each of these choices corresponds to a partition of C objects into P sets 
where order does not count and the objects are indistinguishable. Specifically, for 
a given choice all of the remaining C items are divided among the P sets, where 
everything up to the first of the f — 1 chosen objects goes to the first set, the objects 
between the first and second chosen element goes into the second set, and so on. 
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We may recast the above as saying the number of solutions to xi H \-xp~C, 

with each x, > 0, is ( J_7 ) ■ One of the advantages of this approach is that it is 
very easy to add in lower restrictions. For example, to find the number of solutions 

to yi H \-yp = C with y, > «, for some choice of non-negative integers (with of 

course ni H !-«/>< C), simply let y, = x, + /i,. Now each x, > is integer valued, 

and the number of solutions to the equation in the yis is the number of solutions to 

XI + ...+xp = c-(ni + •••+«/.), or f-("'+;+;''-)+^-'). 

We use the above perspective to analyze the Zeckendorf decomposition of an 
X G [F„,F„+i). The number of such x is F„+\ — F„, which by the Fibonacci recur- 
rence relation is just F„_i. We count the number of x e [F„,Fi,^i) with exactly k+l 
summands in their Zeckendorf decomposition, denoting this quantity by N„{k); we 
choose to record the number of summands as k+l as each x must have at least one 
summand, namely F„. By the standard proof of Zeckendorf s theorem, we know 
each X has at most one such valid decomposition. The existence of a Zeckendorf 
decomposition follows by showing Y!k=o^n{k) = F„_\ (since no number has two 
valid decompositions, it is enough to know that the number of valid decompositions 
equals the number of integers in the interval [f„,/v,+i)). 

Lemma 1.2.2 Let Pn{k) = N„{k)/F„_i, which is the probability an x ^ [PmPn+\) 
has exactly k+l summands iyi its Zeckendorf decomposition. ThenP„{k) = (" j-^ )/^n- 

Proof. If X has exactly ^ + 1 summands, then x = Z^-, + /^, H ^ ^a + ^■t+ 1 ' where 

^i+i =F„,l< ii < (2 < '3 < • • • < ik < 4+1 - and dj := ij - ij^ i > 2 for 2 < y < fc + 1 
(and d\ := /i — > 1). We can recast this in terms of the cookie problem above. 

Clearly di+d2-\ \-dit+i = n (we have a telescoping series and the last summand 

in the decomposition is F„). Let di = xi + 1 and dj = Xj + 2 for 2 < j < k + I. 
Then the number of x that have exactly k+l summands is the number of tuples 

(c/i , . . . , dk+i) with di +d2-\ h d^+i = n, or equivalently the number of tuples 

(xi,. . . ,Xk+i) withxi +X2-\ hx^._|_i =n— {2k + 1). By our combinatorial result. 

Lemma 11.2.11 taking C — n— {2k + 1 ) and P = k+l we see the number of such 
tuples is just 

the lemma now follows. 

Now that we have an explicit formula for N„{k) and Pn{k), all the claims follow 
for the Fibonacci case. We quickly provide a sketch; see IIKKMWI for details. 



Zeckendorf expansion (Theorem [Tli: As each x e [F„,F„^i) has at most one 
Zeckendorf decomposition, the claim follows by counting the number of valid 
Zeckendorf expansions and seeing that this equals F„_ i . This number is just 
Y,"k=Q (" A-^ )■ The summands vanish if A: > [(«— 1)/2J (the binomial coefficients 
are extended so that (^) = if £ > n or £ < 0). We claim 



irr^'-- 
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(1.11) 



We proceed by induction. The base case is clear, and the general case follows 
from using the standard identity that ('") + (^/['j) = ("^/)- Specifically, 



k=Q 



E 

k=0 



n-l- k\ fn-l- k 



" fn-2-{k-\)\ " fn-l-lc 

"^^ fn-2-k\ "fn-l- k 
^ k )+^ 

Fn-2+Fn-l 



(1.12) 



by the inductive assumption; noting F„_2 + F„-i = F„ completes the proof. 

2. Lekkerkerker's Theorem (Theorem|2]i: The claim follows by computing the ex- 
pected value, which is 



l^ik+m,ik) = i+-^j:k^" 



k=0 



«-l *:=0 



(1.13) 



Let 



k=0 



n-l- k 
k 



Straightforward algebra shows 

(g'{n) = {n-2)F„_i-S'{n-2). 
To see this, note 



(1.14) 



(1.15) 
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<g{n) 



n-l-k 
k 



f(„-2-(.-,)) <"-^-<*-')' 

k=l 

L {n-2-e) 



e=o 



^t'l k\{n-\-2k)\ 



{k-\)\{n-2,-2{k-l))\ 
« - 3 - A 



L^J 



(n-2)F„_3-(f(n-2) 



^=0 



«-3-. 



(1.16) 



which proves the claim. Note that we used ( 11.111 ) to replace the sum of binomial 
coefficients with a Fibonacci number 
We study the telescoping sum 



L^J 



£ {^lY{S{n~2(.) + S{n-2{l+\))) 



(1.17) 



f=o 



(which is essentially Sin)). Using (11.15b yields 



L^J 



£ (-l)^'(« - 3 - 2£)F„_3-2£ + 0(F„_2 



(1.18) 



While we could evaluate the last sum exactly, trivially estimating it suffices to 
obtain the main term, which will give us Lekkerkerker's Theorem. 
Binet's formuliU states that F„ = -^- (p" - '^ ■ {\ - (p)'\ with (p = i^ is the 
golden mean. We use this to convert the sum into a weighted geometric series 
(where each factor is multiplied by a simple polynomial): 



S{n) 



-^ £ (n-3-2£)(-l)V-3-2^+0(F„_2) 

I "~^ I I "~^ I 

(«-3) £ (-(p-2)^-2 ££(-(?- Y 



9' 



x/5 



f=0 



1=0 



0(F„_2).(1.19) 



^ It is worth noting that while one can establish Binet's formula by substituting and checking , it can 
also be derived via generating functions; it is no coincidence that in the general proof generating 
functions will play a central role. 
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We use the geometric series formula to evaluate the first term, and differentiating 
identitie^ on the second. After some algebra we obtain 

n(D"^^ n(D" 

As /vi-i — "7? + C'(l) for large n, we finally obtain 

'^{n) = ;^+0(F„_2). (1.21) 

We note a more careful analysis is possible, and such a computation leads to an 
exact form for the mean. 

3. Gaussian Behavior (Theorem O: As we have the density function Pn{k) = 
("^^^ ) /fit- 1 , one way to prove the Gaussian behavior is to show that as « ^- oo 
the function Pn{k) converges to a normal distribution. This may be accomplished 
via Stirling's formula (see IKKMWI for the computation). Another possible ap- 
proach would be to generaUze the proof of Lekkerkerker's theorem to compute 
all moments, and then appeal to the Moment Method; however, as the combi- 
natorics become harder as the moment increases, the Stirling approach is more 
tractable. 



1.3 Positive Linear Recurrence Sequences 

We discuss the main ideas in proving the Gaussian behavior for the generalized 
Zeckendorf decompositions for Positive Linear Recurrence Sequences {//„}. In par- 
ticular, we discuss the obstructions that arise in trying to use the argument from the 
previous section, and describe the techniques that handle them. The computations 
become quite long and technical; we refer the reader to IIMWI for these details, and 
content ourselves with describing the method below. 

The first step is to prove that all integers have a unique, legal representation 
involving the Positive Linear Recurrence Sequence. The proof is essentially just 
careful book-keeping. 

The next step is to determine the size of a general term H„ in our sequence. The 
most important use of this is to count the number of integers we have in a given 
window, which is used to normalize our counts to probabilities. In the Fibonacci 
case, the number of integers in the interval [/s,,^«+i) is just F„+i — F„ = F„_i; we 
then used Binet's formula to approximate well the size of f„_i. This all generalizes 
immediately as we have a linear recurrence relation. 



3 As E"=o^' = ( 1 - -*^"'^' )/{ 1 - ^'). applying -^'s '" 1"°* ^i'l^'^ gi^^^ S"=o ^-'^' =x{\- (in +\)x'" + 
my"+')/(l-x)2. 
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We first sketch another proof of the results from m.2[ and then comment on how 
to generaUze these arguments. While we do not discuss the proof of Theorem|5]in 
detail, it too can be handled by this method (see QMWI ). 



1.3.1 New Approach: Case of Fibonacci Numbers 

We change notation slightly for the rest of this chapter in order to match the notation 
of IIMW] . Let p„ It := #{N <S [F„,F„^i): the Zeckendorf decomposition of A^ has 
exactly k summands}. This double sequence satisfies a nice recurrence relation. For 

any A^ € [F„+i, f„+2)> we have A^ = F„+i + iv H where f < « — 1 (we are not 

allowed to have adjacent Fibonacci numbers in our decomposition). Imagine A^ has 
exactly ^+1 summands in its decomposition. It must have F„+i and it cannot have 
Fn, and needs exactly k more non-adjacent summands from Fi to F„_i. There are 
Pn-i.k ways to have k non-adjacent summands with F„_i included, p„-2.k ways to 
have k non-adjacent summands without F„_i but with F„_2 included, and so on. We 
thus obtain the following formula for p„+i,k+i'- 

Pn+l,k+l ^ Pn-\.k+Pn-2.k-\ ■ (1-22) 

Similarly (replacing « + 1 with n) we find 

Pn,k+\ = Pn-2,k+Pn-3,k-\ ■ (1-23) 

Subtracting, we find 

Pn+ll+l = Pn.k+l+ Pn-l,k- (1-24) 

Our goal is to extract information about the p„ i^. A powerful approach is to use 
generating functions. The generating function in this case is 

L Pn*^'y" = i_/_,„2 - (1-25) 

n.k>0 

Using partial fractions, we find 



1 —y — xy 



l-y-xy^ yi{x)-y2{x) \y-yi{x) y-yiix) 



1 



(1.26) 



where y\ {x) and yiix) are the roots of 1 — y — xy^ = and the coefficient of y" is 
8{^) = lLk>o Pn,kx'' ■ Note the roots are readily computed via the quadratic formula; 
this is not true for the general case, and is the source of much of the technicalities. 

As in the introduction, let K„ be the corresponding random variable associated 
with k. Using the Method of Moments, it suffices to prove the moments of K„ con- 
verge to the moments of a Gaussian to prove Gaussian behavior We do so through 
differentiating identities. With g{x) as above, differentiating once and setting x = I 
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yields 

^(1) = L/Vi = F„+,-F„, (1.27) 

which is just the number of elements in our interval [Fn , /^h+ i ) ■ Differentiating again 
gives 

g'{x) = Y.^pn.,k^-'. (1.28) 

A:>0 

If we take x = 1 we essentially obtain the mean of K„ ; we need to divide by F„_ i as 
the p„,A:'s are counts and not probabilities. As ^(1) = F„_i, we find 

g'{\) = g{l)E[K„]. (1.29) 

We continue, and find 

{xg'{x))' = Y.k^p„,k^-\ (1.30) 

A:>0 

which leads to 

{xg'ix))' = §(!)£[<], (1.31) 



x=l 



and then 



{x{xg'{x))')'l^^ =g{l)ElK'„], 



(1.32) 



and so on. 

Similar results hold for the centralized random variable K'„ = K„ — E[K„]. Miller 
and Wang IMWl prove that £[(/:;,)2'"]/(SD(7i:,^))2'« -> (2m- 1)!! (withSD«) the 
standard deviation of K^;) and E[{K'„f'"-^]/{SD{K;^))^"'-^ -^ 0, which yields the 
Gaussian behavior of A'„. 



1.3.2 New Approach: General Case 

We generalize the arguments from i jl.3.11 to the case of Zeckendorf expansions 
arising from Positive Linear Recurrence Sequences {//„}. We now set /j,,^ = 
#{A^ e [//„,//„+! ): the generalized Zeckendorf decomposition of A^ has exactly k 
summands}. We can find a recurrence relation as before. For the Fibonacci num- 
bers, we found 

Pn+l.k+l ~ Pn.k+l+ Pn,k- (1.33) 

In the general case, we have 

Pn+l.k = Y, H Pn-m.k-j, (1.34) 

m=0 ,;=.v,„ 

where ^o = 0, s,„ = ci + C2 H h c,„ (with these c's the same c's as in the definition 

of our recurrence sequence). 
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The generating function in the Fibonacci case was relatively straightforward, be- 
ing ^ ^ j . Now our generating function equals 

We again perform a partial fraction expansion in order to glean information about 
the coefficients. Instead of 



y / 1 1 





yiix)-y2ix) ' 


\y--yi{x) y-yiix))' 


now we have 








1 ^ 


B{x,y) 
~yiix))nj^iiyj{x)-yi{x)r 


where 

B{x 


n<L 


L-l.«;„+l-l 

m=0 j=Stii n<L—m 



(1.36) 



(1.37) 



(1.38) 

and the y,(jt:)'s are the roots of A(3;) = I -'L^nJoL'f^.'^^'' xh'"""^ = 0. One of the 
major difficulties of the proof is the analysis of the roots yi{x). Unlike the Fibonacci 
case, where it was easy to write down simple expressions for these via the quadratic 
formula, in general the arguments become quite involved. If however the coefficients 
in the recurrence relation are non-increasing then the proofs are easy; see Appendix 
C of IMWI . The general cases (see Appendix A of IMWI ) is more complicated, 
involving continuity and the range of the |3',(x)|'s. The main idea is to first show that 
there exists x > such that A (y) has no multiple roots and then prove that there are 
only finitely many x > such that A (y) has multiple roots. 

The coefficient of y" is g{x) = Y.kyoPn.kX^- We use the method of differentiating 
identities as before to get the moments, and find that K„ converges to a Gaussian as 



The proof of Theorem |5] proceeds similarly. We state the two key elements, the 
recurrence relation and the generating function, without proof, and refer the reader 
to OMWI for the details. Let Pn.k.e be the number of far-difference representations 
of integers in {S„-\,S„] with k positive summands and £ negative summands. The 
recurrence relation is 

Pn,kl = Pn-l,kI+Pn-4,k-lI + Pn-3,e,k-U n>5, (1-39) 

and the generating function is 
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g{x,y,z) = Y. /'«,mAV. (1.40) 

n>0,A:>0,^>0 



1.4 Conclusion and Future Research 

The combinatorial approach has extended previous work, allowing us to prove Gaus- 
sian behavior for the number of summands for a large class of expansions in terms 
of solutions to linear recurrence relations. This is just the first of many questions 
one can ask. Others, which we hope to return to at a later date, include: 

1 . Are there similar results for linearly recursive sequences with arbitrary integer 
coefficients (i.e., negative coefficients are allowed in the defining relation)? 

2. Lekkerkerker's theorem, and the Gaussian extension, are for the behavior in in- 
tervals [F„,F„+i). Do the limits exist if we consider other intervals, say [F„ + 
gi{F„),F„ + g2{Fn)) for some functions gi and g2? If yes, what must be true 
about the growth rates of gi and ^2? 

3. For the generalized recurrence relations, what happens if instead of looking at 
^"^jfl, we study ^"^jmin(l,a,)? In other words, we only care about how many 
distinct //,'s occur in the decomposition. 

4. What can we say about the distribution of the largest gap between summands in 
the Zeckendorf decomposition? Appropriately normalized, how does the distri- 
bution of gaps between the summands behave? What is the distribution of the 
largest gap? How often is there a gap of 2? 
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